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ABSTRACT 


Stationary  stable  processes  that  are  Fourier  transforms  of 

symmetric  stable  independent  increments  processes  are  shown  to  have 

a.s.  finite  conditional  expectation  of  Xt  given  Xg  and  conditional 

variance  of  X  given  X  c.  X  The  associated  conditional 

t  t-o  t-2o 

expectation  predictors  are  nonlinear  in  {Xg,  s<t}  but  are  mixtures 
of  predictors  of  the  usual  type  based  on  the  Gaussian  model. 
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I  INTRODUCTION 


The  work  presented  here  is  part  of  a  program  extending  the 
classical  theories  of  prediction,  detection,  and  smoothing  of 
signals  to  encompass  models  in  which  Gaussian  noise  is  replaced  by 
symmetric  a-stable  (SaS)  noise,  where  0  <  a  (  2  is  the  index  of 
stability.  This  is  consistent  with  the  full  implications  of  the 
theory  of  errors  (Levy,  1925)  and  includes  Gaussian  noise  as  the 
case  a  =  2. 

It  is  important  to  realise  that  Gaussian  noise  (a  =  2)  results 
in  the  frequency  domain  models  being  identical  with  the  time  domain 
models,  by  Bochner's  Theorem  (Doob.  1953).  But  the  same  is  not  true 
for  a  <  2.  That  is,  with  the  notable  exception  of  a  =  2,  the  class 
FT-SaS  of  processes  which  are  Fourier  transforms  of  SaS  noise  (e.g. 
wave  motions)  will  not  be  the  same  as  the  class  of  processes  that 
are  the  outputs  of  linear  systems,  such  as  ARMA  models,  driven  by 
SaS  noise  (Cline  and  Brockwell,  1985),  (Makagon  and  Mandrekar, 

1987). 

Substantial  progress  has  been  made  on  the  prediction  problem 
for  processes  of  the  FT-SaS  type,  i.e.  stationary  processes  = 

/  e**"^  Z(dX)  where  Z  is  an  independent  increments  SaS  process.  In 


the  course  of  this  work  a  number  of  basic  methods  have  been 


developed  which  also  have  application  to  processes  of  the  latter 
type,  including  ARMA  models. 

To  understand  the  basic  approach  of  this  paper  it  is  important 

to  know  that  2n  i.i.d.  SaS  random  vectors  tend  to  separate  out  like 
—  1  //> 

+  j  .  1  <  j  i  n  (LePage,  Woodroofe,  and  Zinn,  1981).  and  are  in 
general  representable  as  mixtures  of  Gaussian  distributions  with 
differing  covariances  (LePage,  1980,  appended  to  this  report).  This 
allows  us  to  think  of  i.i.d.  stable  noise  as  inhomogenous  Gaussian 
noise  with  random  a/2  stable  covariances.  As  might  be  inferred  from 
remarks  made  above,  the  consequences  of  this  inhomogeneity  when 
noise  is  in  the  form  of  wave  amplitudes  differ  greatly  from  the 
effects  of  noise  entering  in  the  form  of  terms  driving  difference 
equations. 

In  recent  work  (Makagon  and  Mandrekar.  1987)  there  is  defined 
the  concept  of  a  generalized  spectrum  for  any  strictly  stationary 
stable  process.  In  case  the  spectrum  is  given  by  independently 
scattered  measure,  as  above,  they  have  shown  that  one  can  obtain  the 
1 Inear  analysis  of  such  signals  by  methods  analagous  to  those  used 
for  Gaussian  processes.  This  extends  (Cambanis  and  Soltani,  1982) 
to  the  case  a  i  1  by  a  more  general  method. 

However.  1  Inear  prediction  is  not  in  general  as  good  as 
conditional  expectation  prediction.  The  latter  has  not  been  much 


studied  for  stable  processes  because  of  the  fact  that  E  | X t |  is 

infinite  for  a  <  1.  Surprisingly,  by  exploiting  (LePage,  1980)  we 

will  prove  for  FT-SaS  processes  (Xt)  that  for  every  a  the  condi tonal 

expectations  E( |Xf |  |  Xg)  and  E(X£  |  X£  g,  X g)  are  almost  surely 

f ini te.  This  is  surprising  in  view  of  the  fact  that  for  a  <  1  the 

expectations  E  | X ^ |  are  infinite,  and  for  all  a  <  2  the  expectations 
2 

E  X£  are  inf  ini te . 

As  a  consequence  of  these  conditional  moments  existing,  the 
conditional  expectation  predictor  of  X^.  regardless  of  the  number  of 
predicting  variables,  is  well-defined  and  optimal  for  conditional 
mean  squared  error,  provided  the  predicting  variables  include  two 
time  points  of  the  form  t-6,  t-26  for  some  6^0.  By  combining  these 
results  with  (Cambanis  and  LePage  1987)  it  can  be  shown  that  for 
a  <  2  the  conditional  expectation  predictor  EfX^  |  Xt_g,  <5q<6<L)  is 
asymptotically  consistent  for  X(  as  6  -*  0  (excluding  0  in  the 
discrete  case)  and  L  -»  “,  whereas  the  linear  predictor  is  not. 

Since  the  above  conditional  expectations  turn  out  to  be 
a-posteriori  averages  of  Gaussian  condi tonal  expectations  computed 
for  various  covariances,  these  non-linear  conditional  expectation 
predictors  are  a  smoothing  of  Hilbert  space  methods  and  are  in  fact 
Bayesian  predictors  for  a  naturally  occurring  a-priori  distribution 


intrinsic  to  FT-SaS  processes. 


The  new  methods  generally  allow  us  to  compute  quantities 
previously  thought  to  be  undefined  because  stable  r.v.  lack  certain 
moments.  For  example,  the  identity  |  Xg)  =  E  cos(A(t-s)), 

where  A  is  a  random  sample  from  the  normalized  spectral  distribution 
function  was  known  for  1  <  a  £  2  {Kanter  and  Steiger,  1974). 
Corollary  2.1.2  below  proves  this  result  for  0  <  a  <  2  by  a  totally 
new  direct  calculation  which  does  not  require  existence  of  the 
unconditional  expectation  and  bypasses  differentiation  of  the 
characteristic  function  altogether.  Certain  other  conditional  and 
unconditional  integrals  can  be  directly  calculated  by  the  same 
method,  including  the  integer  moments  of  the  characteristic  function 
of  the  processs  conditioned  on  the  invariant  sigma  algebra  (Cambanis 
and  LePage,  1987). 

As  mentioned  previously,  this  work  is  based  on  a  representation 

of  FT-SctS  processes  as  mixtures  of  stationary  Gaussian  processes 

(with  randomly  chosen  covariance  function  9)  due  to  (LePage,  1980). 

The  new  observation,  specialized  to  the  case  3  =  o{X Xq)  and  X  = 

Xj.  is  that  the  conditional  density  of  (X^ ,  Xq) ,  given  the  sigma 

field  generated  by  the  covariance  function  9,  cancels  terms  in  the 

conditional  expected  squared  error  of  prediction.  This  forces 

convergence  of  the  conditional  expectation  E?(X-E?X)2.  Such 

3 

integrals,  including  E  X,  are  computed  as  mixtures,  on  9,  of 


6 


;jr«.rw 


9-specific  Gaussian  integrals.  The  indicated  predictors  are  thus 
conditional  mixtures  of  0-specific  Gaussian  linear  predictors. 

STABLE  PROCESSES 

The  log  characteristic  function  of  any  stationary  Gaussian 

random  functions  having  a  continuous  covariance  takes  the  form 

log  E  exp  i  ^=1rkXt  =  -(o2/ 2)  El^  ^  exp  it^J2.  (1.1) 

k 

2 

where  a  >0  and  is  a  random  variable  whose  probability 
distribution  function  is  the  normalized  spectral  distribution 
function.  Stable  analogues  of  these  laws  may  be  obtained  by 
replacing  the  exponent  2  in  (1.1)  by  a  number  a  in  the  range  0<a<2. 
In  (LePage,  1980)  it  was  proved  that  the  resulting  characteristic 
functions  are  precisely  those  of  the  class  of  FT-SaS  processes.  The 
following  construction  for  such  processes  was  given.  For  each  0<a<2 


define  r.v.  (X(t).  telR)  by 


X  =  2.  lCos(A.t  +  ®.)Y.rV 
t  J=1  v  J  J  J  j 


(1.2) 


In  (1.2).  which  converges  a.s.  for  each  t,  the  sequence  of  r.v.  { A^. } 


(which  we  denote  by  A  for  brevity)  are  i.i.d.  from  any  distribution 


on  IR;  ®  are  i.i.d.  uniformly  distributed  on  the  interval  [-7r,rr]:  Y 
are  i.i.d.  with  E|Yj|a  finite;  T  are  the  consecutive  arrival  times 
of  a  homogeneous  Poisson  process  with  unit  intensity  function  on  the 
time  domain  IR+;  and  the  sequences  of  r.v.  A,  ®,  Y,  f  are  mutually 


independent . 
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From  (LePage,  1980,  Theorem  7.4)  the  series  (1.2)  converges 

almost  surely  for  each  telR  and  0<a<2,  and  the  r.v.  defined  by  the 

left  side  of  (1.2)  have  log  characteristic  function 

2_aB(a)C(a)  EIyJ*  E|s£=1  rk  exp  itjAj®.  (13) 

where  B(a)  =  a  J*  (cos(r)-l)  — and  C(a)  =  |l+e^|adTj. 

r 

From  (1.3)  it  may  be  seen  that  the  law  of  X  depends  on  the  law 

of  Y.  only  through  the  a-th  absolute  moment.  Taking  (Y.)  to  be 
•  J 

Rayleigh  distributed,  equivalent  to  letting  { Y ^ }  be  complex  Gaussian 
and  taking  the  real  part  of  (X(t),  telR) ,  yields  a  process  X  which  is 
conditionally  Gaussian  and  stationary  given  the  sequences  A,  T.  The 
process  given  by  the  infinite  series  (1.1)  may  be  written 

Xt  =  JR  eiXt  Z(dX).  (1.4) 

where  Z  is,  conditionally  on  A,  T,  a  Gaussian  orthogonal  random  set 
function  supported  on  the  sequence  of  frequencies  A  and  -A,  with 

Z((-Ak>)  =  Z((Ak})  and  Z((Ak>)  given  by 

=  (|)  2j=1rj1/aY.(I(Aj=Ak)ei®j  +  I(-Aj=Ak)e_i®j),  k*l.  (1.5) 
From  (1.2)  or  (1.5)  it  is  seen  that  the  conditional  covariance 
function  of  X  given  A,  T,  is  given  by 

qo  -O  / 

0(t)  =  2  irj  cos(Ajt)/2  telR.  (16) 

That  is,  (X(t),  telR  is  conditionally  Gaussian  and  stationary  with 


covariance  (1.6)  given  A.  f. 


2.  EXISTENCE  OF  CONDITIONAL  MOMENTS 
gj  i  3 

We  require  the  formula  E  rj  =  JV(d0p)  EgT),  where  r?  is  a 

non-negative  polynomial  in  the  coordinate  r.v.  of  the  random 

function  X,  3  is  the  o-field  generated  by  a  finite  number  of  the 
3 

coordinate  r.v.,  Egi7  is  the  Eg-condi tional  expectation  given  3,  and 

Tr(d0|?)  is  the  regular  conditional  probability  distribution  of  the 

3 

finite  number  of  coordinate  values  of  0  which  appear  in  EqH,  A 
proof  may  be  given  by  martingale  methods  from  the  corresponding 
result  for  the  discrete  case  which  is  P(A|D)  =  20  P(0AD)/P(D)  =  20 
P(0|D)  P0(A|D). 

Theorem  2.1.  Let  X  be  a  random  function  which  is  a  mixture  of 
zero-mean  Gaussian  stationary  processes.  For  real  numbers  s,t  the 
conditional  expectation  E  | X 1 1  is  almost  surely  finite,  where  3  is  a 
CT-algebra  with  respect  to  which  X^  is  measurable. 

x0 

Proof .  It  is  enough  to  prove  that  E  | X ^ |  is  finite  a.s.  If  0(0)  = 

0,  then  XQ  =  Xj  =  0  almost  surely.  For  each  fixed  0  with  0(O)>O, 

the  conditional  distribution  of  X^  given  Xq  is  normal  with  mean 

X0  2 
4(0,0, 1)  =  E_  X.  =  (0(l)/0(O))Xn,  and  variance  o  (0,0,1)  = 
y  i  u 

0(O)(1-(0(1)/0(O))2).  Then. 

X  X 

E  |X(1)|  -  -T^0^OpOj  E0  |Xj|  Tr(d0|xo) 

S  -T(e(o)>0)  (iM(B.O.l)l  ♦  *(0,0,1))  *(d0|Xo).  (2.2) 


Since  1 0 ( 1 )  |  ^ 1 0 (0) |  a.s.,  the  p-term  of  (2.2)  is  a.s.  finitely 

integrable.  Now  use  rr(d0  I^q)  =  ir(d0)Tr(XQ |0)tt(Xq)  .  where  ir(r|0)  = 
-1/9  -1/2  9 

(2ir)  0  (0)  exp(-r  /20{O)),  from  which  it  is  seen  that  the 

cr-term  of  (2.2)  is  also  a.s.  finitely  integrable.  □ 

Corollary  2.1.1.  If  ?  is  a  sigma  algebra  generated  by  a  subset  of 
the  random  variables  (X(t),  te!R} ,  then  E^X^  =  J*  Tr(d0|l?)  E^X^.  which 
is  a  mixture  of  the  0-specific  Gaussian  best  linear  predictors  of  X^ 

X0 

based  on  9.  In  particular,  E  Xj  =  E(0( 1 )/0(O) jXQ)  XQ  =  (X(8(0j>0) 
ir(d0|Xo)(0(l)/0(O))  XQ.  □ 

Corol lary  2.1.2.  (Cambanis-LePage) .  For  a  process  of  the  form  (1.2) 
with  0<a£2,  E(0(1)/0(O) |XQ)  =  E  cos(A1). 

Proof .  The  case  a=2  is  obvious.  For  a  <  2,  from  (1.2)  and  (1.6) 

x.  =  j"  cos(0.)  rT1/2 
0  J=1  j'  J 

0(1  )/©(o)  =  2°.  r~2/a  cos(x . )/2°?  .  r~2/a. 
j=i  j  r  j=i 

Since  A,  0,  T  are  independent  and  Xq  is  measurable  a{G,r}. 

X  X 

E  °  (0(1)/0(O))  =  E  0  E(0( 1 )/0(O) | 0 , r ) 

Vs  r-2/a  „  , .  >.  “  r-2/a 

=  Sj_,  Tj  E  costAj)/!^,  Tj 

=  E  ros( A^ ) .  □ 

Since  no  special  role  is  played  by  times  zero  and  one,  the  formula 
E[0( t)/0(O) |X(s)]  =  E  cos(A^(t-s))  follows  at  once.  This  result 
extends  (Kanter  and  Steiger,  1974)  to  the  case  a  £  1. 
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Proposition  1.  There  exist  FT-SaS  processes  for  which  the 


0  0  2 

conditional  mean  squared  error  E  (X^-E  Xj)  is  infinite. 

Xq 

Proof.  Let  H(Xq)  =  J  (0( 1)/0(O))  7r(d0|Xo).  That  is,  E  u  X]_  =  H(XQ) 

Xq .  Then 

X  X 

E  °(XrE  °X1)2  =  X(e(0)>0)  6(0)  (1-(0(1)/0(O))2)  7r(d0|Xo) 

+  X(0(O)>O)  X0  (H(XO)-(0(1)/6(O)))2  ir(d0 |XQ) .  (2.3) 

The  second  term  on  the  right  of  (2.3)  is  integrable.  The  first  term 
on  the  right  may  be  handled  using 
ir(d0  |XQ)  =  ir(X0(6)  7r(d0)/ir(Xo) 

-1/2  -XO/20(°) 

=  (2tt0(O) )  1/Z  e  U  7r(d0)/ir(Xo).  (2.4) 

From  (1.6)  it  is  seen  that  0(0)  possesses  a  stable  distribution  of 

index  a/2  under  tr(d0).  Taking  with  support  in  the  two  point  set 

{-tt/2,  tt/2}  ensures  0(1)=O  a.s.,  in  which  case  the  integrand  of 

(2.3)  is,  by  (2.4),  asymptotically  of  order  r  at  infinity. 

This  fails  to  be  integrable  if  a£l.  □ 

On  the  other  hand,  the  conditional  expected  squared  error  of 

prediction  given  3  =  o{\ ,  Xq)  is  a.s.  finite. 

Theorem  2.2.  Let  X  be  a  random  function  of  the  form  (1.2).  For 

real  numbers  r.s.t  with  |r-s|  =  |  t-s |  the  conditional  expectation 
5  5  2 

E  (X^-E  Xt)  is  almost  surely  finite  if  9  is  a  o-algebra  with 
respect  to  which  X^  and  Xg  are  measurable. 


j  win  w  m  mm  w  w  m  w  w  u  W  WWWU  m  l  JJ 1 1  f  1 1 1 1  n  1 1 )  i )  u  i.  gygy  L«  w  n^ynr 


Proof.  Replace  7r(d0  |XQ)  by  Tr(d0  |X_1  ,XQ)  =  Tr(d0|5)  in  (2.3),  and 
denote  by  p  the  correlation  p  =  0(1)/0(O).  Then, 


T(d0l?)  ^  (ir(d0)/ir(X_1 ,  X0))(2t r)  10~1  (0) ( 1-p2)"1 
From  theorem  2.1. 

E?X2(1)  =  X  7T(d0|?)  EgX2 

=  s  ir(d0|5)  [E^(X1-E^X1)2  +  (EgX1  )2] 


(2.5) 


(2.6) 


^ie  first  term  in  the  brackets  of  (2.6)  is  less  than  or  equal  to 

0(0)  (1-p  ),  and  is  therefore  ir(d0)-integrable  when  multiplied  by 

the  Jacobian  from  (2.5).  The  second  term  in  the  brackets  of  (2.6) 
2 

may,  on  (p  ^1)  be  written 

[P  X  +  -i(.2)  -  P0^1)  (X  -  p  X  )]2.  (2.7) 

0(0)  (1 V)  10  v  ; 

The  first  component  in  (2.7)  is  p  Xq  whose  square  is  integrable  by 
Tr((d0|?).  It  is  therefore  enough  to  bound 

D(0)  .  pta  -  Pfl(U  |  (2.8) 

1  6(0)  (1-p2)  J 

^■et  II  II  denote  Lg  norm  with  respect  to  the  zero-mean  stationary 

Gaussian  process  with  covariance  function  0.  Then 
9(0)  =  I  IXj  1 12  >  |  |pXQ|  |2  +  |  | D( 0 )  (X_rpXQ)  I  |2 
=  P20(O)  +  D2{0)  0(0)  (1-p2). 

Therefore  D2(0)  <  1  on  p2/l .  But  on  p2=l  (2.7)  is  replaced  by  XQ. 


Thus  (2.5)  is  a.s.  finite.  Note  that  the  Cauchy-Schwar tz  bound  on 
(2.8)  is  inadequate.  □ 

3.  FUTURE  DIRECTIONS 
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The  ergodic  properties  of  {Xf}  are  now  completely  understood. 

In  (Cambanis  and  LePage,  1987),  subject  to  mild  conditions  on  the 
spectral  distribution,  a  complete  characterization  of  the  invariant 
sigma  field  of  the  process  {Xt}  is  obtained  in  terms  of  amplitude 
and  frequency  variables  measurable  with  respect  to  the  remote  past. 
The  non-ergodic  component  is  identified  with  i.i.d.  uniformly 
distributed  phase  variables. 

Using  the  above  results,  it  follows  that  under  rather  weak 
conditions  the  predictors  E(Xt  (  X  g,  6q<6<L)  converge  almost 
surely  to  X^  (and  condi tional ly  in  mean  square  given  the  random 
spectral  amplitudes)  as  L  -»  °°,  6q-»  0.  Thus,  the  action  of  SaS  noise 
in  wave  amplitudes  is  to  make  prediction  ultimately  more  perfect 
than  would  be  the  case  for  Gaussian  noise.  Since  by  a  result  of 
(Makagon  and  Mandrekar,  1987)  linear  prediction  optimized  for  the 
given  FT  SaS  distribution  is  only  consistent  for  Xt  under  very 
exceptional  conditions  on  the  spectral  measure,  it  follows  that 
linear  prediction  is  inferior  to  conditional  expectation  prediction 
for  a  <  2.  Such  results  generalize  easily  to  the  multiparameter  and 
dim  >1  cases,  and  will  apply  to  certain  types  of  spatial  processes 
which  are  Fourier  transforms. 

By  contrast  with  the  above,  linear  prediction  for  linear  system 
SaS  processes,  such  as  ARMA(p.q)  models  driven  by  i.i.d.  SaS  random 
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vectors,  is  intrinsically  interesting  due  to  the  way  in  which  the 
noise  enters  additively.  The  conditional  expectations  E( |X^ j  |  , 

s  <  t)  are  a.s.  infinite  except  in  trivial  cases  due  to  E( |e^ |  |  X  . 
s  <  t)  =  «,  i.e.  the  conditional  expectation  of  the  independent 
error.  By  Theorem  2.1  it  follows  that  such  processes  cannot  be  of 
the  FT  SoS  type. 

If  it  can  be  proved  that  E ( | J  |  X  ,  s  <  t)  <  “  a.s.  for  all  r 
<  t,  then  an  obvious  choice  for  non-linear  prediction  will  be  to 
predict  the  value  zero  for  €  and  use  the  conditional  expectation 
predictor  on  r  <  t.  This  most  promising  approach  is  currently 
under  study. 


/.  r y,  ^y.  * ,  * ,  • . « 4  y,» .  • ,  *\  *  .  *.  •*.  /.  -• 
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MULTIDIMENSIONAL  INFINITELY  DIVISIBLE  VARIABLES 
AND  PROCESSES.  PART  I:  STABLE  CASE 


By 

Raoul  LePage 


1 .  Introduct ion. 

Series  deconposit ions ,  involving  the  arrival  tines  of  a  Poisson 
process,  have  been  given  by  Ferguson  and  Klass  [1]  for  the  non-Gaussian 
component  of  an  arbitrary  (real-valued)  independent  increments  random 
function  on  the  unit  interval.  LePage  and  Uoodroofe  and  Zinn  [4]  have 
rediscovered  a  variant  of  this  decomposition  in  connection  with  their 
study,  via  order  statistics,  of  the  limit  distribution  for  self-normalized 
sums  (e.g.  Students  ’ -t ) ,  when  sampling  from  a  distribution  in  the  domain 
of  attraction  of  ar.  arbitrary  stable  law  of  index  •'  2. 

The  present  paper  obtains  a  cha ratter izat ion  of  static  laws  cn 
spaces  of  dimension  greater  than  one.  This  characterization  is  formally 
like  that  of  Ferguson-Klass  for  dimension  one,  but  with  i.i.d.  vector 
multipliers  or.  the  Poisson  terms.  The  law  of  .these  coefficients  may  be 
chosen  proportional  to  the  Levy  measure,  although  this  is  not  necessarv. 
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function,  and  (iv)  how  to  construct  an  arbitrary  harmonizable  stationary 
symmetric  stable  random  function  having  multidimensional  domain  and/or 
range.  Symmetric  stable  laws  are  shown  to  be  mixtures  of  Gaussian  laws. 

Partly  because  of  the  self-contained  character  of  Kuelbs'  paper 
[3],  in  which  the  characterization  of  the  log-characteristic  function 
of  a  stable  law  is  extended  to  real  separable  Hilbert  space,  the 
Hilbert  space  level  of  generality  has  been  chosen  for  this  paper. 

Later  extensions  of  Kuelbs'  result  to  3anach  and  more  general  spaces 
support  a  corresponding  generalization  of  these  results.  In  addition 
to  Kuelbs'  result  we  need  a  method  employed  by  Ferguson  and  Klass  to 
transform  certain  dependent  series  into  eventually  identical  indepen¬ 
dent  ones.  We  also  require  standard  results  giving  conditions  under 
which  an  independent  series  in  Hilbert  space  converges  almost  surely 
(e.g.  [2],  Theorem  5.3).  The  rest  of  the  paper  is  basically  self- 
contained  and  affords  a  surprisingly  accessible  and  clear  view  of  :  <  2 
stable  laws,  and  random  functions,  based  on  elementary  series  constructions. 

Part  II  of  this  paper  will  generalize  these  results  to  the  infinitely 
divisible  case. 


2. 


Notation. 


The  following  symbols  and  conventions  will  be  in  force  throughout 
this  paper. 


(2.1) 


n. 

A 

ac 

D 

3 


=i> 

a 


H 


"is  asymptotic  with” 

"equals  by  definition" 

"has  the  same  distribution  as" 

"converges  in  distribution  to" 

0  <  a  <  2 ,  an  index  of  stability 

arrival  times  of  a  Poisson  process 
with  unit  rate 

a  real  separable  Hilbert  space 


The  material  of  the  next  section  is  drawn  from 
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3. 


Motivation. 


Limit  theorems  are  not  *he  subject  of  this  paper.  However,  we 
should  not  proceed  without  benefit  of  the  following  example,  which 
exposes  some  connections  between  a  <  2  stable  r.v.  and  the  Poisson 


process. 

Let  {e., j  >  1}  be  independent  of  the  sequence  T  and  i.i.d. 

with  PCe^  =  1)  *  P(el  *  -1)  *  j,  and  define  G(x)  =  x-01,  Vx  1  . 

Think  of  G(x)  “  P(jx:  >_  x)  Vx  1,  where  X  is  a  r.v.  symmetrically 

distributed  about  zero.  We  will  construct  particular  r.v.  X^,...,X  , 

i.i.d.  as  X,  whose  normalized  suns  converge  in  distribution  to  the 

symmetric  stable  law  of  index  a.  To  do  this,  use  the  arrival  times  of  a 

__  1 

Poisson  process  to  generate  uniform  order  statistics,  apply  G  to  these, 
multiply  by  the  signs  t ,  and  permute.  As  constructed,  the  normalized 


«-  ^  —l/U 

sums  will  actually  converge  almost  surely  to  C  :  ^  (see  (3.1) 


below) ,  a  series  possessing  the  symmetric  stable  law  of  index  n. 

~ao  —  1  /  cx 

In  fact,  a  direct  uroof  of  the  stability  of  1.  E.f.  follows 

1  1  J 

easilv  from  the  observation  that  the  arrival  times  of  several  (say 


K  >  1)  independent  unit  rate  Poisson  processes  (run  simultaneously) 
constitute  K_1  times  the  arrival  times  of  a  unit  rate  Poisson  process. 
This  argument  works  just  as  well  for  £  replaced  by  any  vector  sequence 
(provided  the  series  converges)  and  suggests  the  multivariate  extensions 
of  sections  4  and  5. 


For  each  n  >  1  let  U.,.  <  l'.  <  •••  <  U,  .  denote  the 

(1)  -  (2)  -  -  (n) 

order  statistics  of  i.i.d.  random  variables  U. , . . . ,U  which  are 

I  n 
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uniformly  distributed  on  [0,1].  Then  for  each  fixed  n  >  1,  letting 


X.  =  e.G  1(U.),  Vj  >  1, 
1  1  J  ~ 


(3.1) 


-1/a  _n  „  D  -1/a  _n  -1  . 

n  X,  e.X.  =  n  _  e.G  (U,..) 
1  J  3  11  (l) 


-1/a  vn  TT-l/a 

n  X  £.U 

1  1  (l) 


D  -1/a  „n  .  .  v-l/a 

-iWW 


-  i"  E,r.'1/:,(r 

1  33  n+1 


1.  j 

4a.  s.  -i-a . s .  (SLLN) 

eT  t.rTVa  i 

1  1  1 


It  is  convenient  to  refer  to  e^Ih,  J  1  ^  as  c^e  residual  order  statistics, 

keeping  in  mind  that  the  ordering  is  on  decreasing  absolute  values. 

The  same  example  suggests  an  invariance  principle  (proved  in  [4]) 

for  self-normed  sums  such  as  X1?  X.  /fcx2  which,  regardl  ess  of  a  , 

1  j  1  n 

converge  in  distribution  to  a  limit  law  depending  only  on  the  stable 
attracting  X^.  For  the  r.v.  constructed  above,  Vn  1, 


1.  •  oo  _2  /  ct 

Use  I\  n,  j  and  X^  j  <  00  a.s.,  and  apply  the  3-series  theorem 

conditional  on  the  sequence  f. 


I_  «  +  m  «_*'  .  •  «  /  •  *  •  '  •  '  »  *  *  *  * 


-  •  N  «  ■  »  . 

'/  .*  V,  /  /,  .  . 

•.  *.  .  I  ’  • 


p 


;«  *,■  *•*_,»* 
*  ■  t  *  .  *  •  .  “  - 


(3.2) 


d 


1  -1  .1 _  a-s- 


00  -l/: 

l  e  r 
1  i  j _ 


lzn  x2 

/z"  r?2/:i 

k:  t:2^ 

J  1  j 

<  i  i 

J  1  3 

2 .  . 

That  is,  the  limit  law  of  the  t-statistic  is  that  of  the  t-stat istic 

calculated  on  the  residual  order  statistics  (see  also  [6]). 

Even  the  construction  of  stable  independent-increments  processes 

can  be  motivated  by  means  of  the  same  example.  We  restrict  our  atten¬ 


tion  to  the  homogeneous  increments  case.  Let  [T  ,j  1  ^  1. 


and  suppose  the  sequences  T,  e,  T  are  mutually  independent.  The 

partial  sum  processes  r[nt]  X  ,  0  <  t  <  1,  Vn  >  1  can  be  effected  by 
independent  selections  of  X^,...,X^  into  subsets  of  sizes  [nt]  using 
multiplication  by  indicators: 


(3.3) 


I^n)(t)  =  I(T1  ± 


.  ,  .  [nt]  - 1  .  IU)(t) 

i-  (t)  a  kt.  < - ^ri^T - } 


1  £  j  <  n 


Then  for  each  n  >  1, 


(3.4) 


{n 


-1/a  v[nt] 


r{nt,x,,  tcto.il!  S  !Zj  tc[0-11 

v  a. s.  in  D[0,1] 


.  1/a 


{i*  i (t .  <  t)-.:.1/a,  tc[o,i] / 

i  j  -  jj 


Details  of  this  argument  are  unpublished. 


2*  The  square  of  this  t-statistic  is  simply  related  to  the  square  of 
Students’-t,  and  both  have  the  same  limit  law. 
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V, 

V, 


>5 


>3 


hi 


,y 


■n 


4. 


Stable  Laws  on  H. 


Suppose  x , j  >_  1;  are  i.i.d.  random  vectors  in  a  real  separable 

Hilbert  space  H  and  that  the  sequences  X,  T  are  independent.  For 

each  n  >  1  denote  by  K  the  number  of  arrival  tines  T  in  the 
—  n 

^  xx  JL 

interval  [0,a  ],  a  =  I.  j  .  This  choice  of  a  is  from  [11.  Its 
n  n  x  xx 

advantages  will  be  apparent  in  what  follows. 

K 

Remark.  Sums  of  the  kind  E  n  are  for  each  n  >  1  defined  to  zero  on 

the  event  K  =  0.  Use  (  ,  ),  |J  ||  ,  to  denote  H  inner  product  and  nom 
n 

For  each  n  >.  1,  xcH,  c  >  0,  (see  also  [1],  pg .  1639), 

i(x, cl  nX  ) 

(4.1)  Ee  1  j  J 


Kuelbs  ([3],  lemma  2.2)  has  proved  that  the  log-characteristic 
function  of  a  (non-Gaussian)  stable  probability  measure  on  H  is 
necessarily  of  the  following  form,  for  a  unique  S c  H  and  finite 
Borel  measure  0  on  S  =  { xeH I ! ■  x =  1; , 


k  id*  n(x,x.)r:1/;t 

EE  ne  1  ^  1 


,  r  a  ic (x ,X  ) t  K 

E(a-1  n  Ee  1  dt)  n 

n  J0 


fa.  ic(x,X  )t“1/-i 

E  J0  (e  1  -l)dt 


a  I  w  ,  i(x»Xl)r1N  dr 

Ec  a  /  -1/a  (e  -1) 

n  r 


(4.2) 


-  1  - 


i(x,g) 


+  f  r  (ei(x’s)r 

JsJo 


i(x, s)rv  dr 


) 


1+r 


1+a 


a (ds)  ,  Vx  z  h 


It  is  convenient  to  refer  to: 

(4.3)  expression  (4.2)  with  6=0. 
Define  6  =  g( S) ,  y  =  a/5,  and 


0 


n 


dr 

1-KX 

r 


a (ds) ,  Vn  >  1  . 


Lemma  4 . 4  If  (4.3)  is  the  characteristic  function  of  a  stable  law  on. 
H  and  if  the  sequence  (X^  ,j  j>  l)  is  i.i.d.  y  and  independent  of  the 
sequence  T  ,  define 

x(n)  A  a-l/afil/“  z  n  x  Vn  >  1  . 

1  J  J  n 

Then  VxcH,  E  exp  Kx.X^)  converges  to  (4.3). 

Remark .  This  result  is  not  altogether  satisfactory  since  convergence  of 
the  series  {x^n\  n  >  1}  is  through  stochastic  times  {K  ,  n  >  1;  , 
and  is  not  yet  a.s.  in  H.  These  defects  are  remedied  in  Theorem  4.8 
below. 

Proof .  For  each  x  e  H,  n  _>  1,  by  (4.1), 
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t.  ...  v  v  a.-v v  v  v ■■■ vv -v  v viv/rS-S-;. 


(4.5)  log  Eel(x,X  }  =  - 


i(x,8  ) 
n 


f°°  i(x,X  )r 

i(,-1/^1/aa-1/ci(e  _1) 


•  IU 


'a  <s  '  a 


/-i(x,s)r  ,  i(x,s)rv  dr  N 

1/ct  -l/a(e  -i-— c(ds) 

an  1+r  r 


(4.3)  q 


From  ([1],  lemma  2),  we  conclude  that 


f<-Kn+l  r-l/i  .. 

*K  Vj  • n  - 1 


are  independent.  Furthermore,  using  an  argument  drawn  from  [4],  Vn  >  1, 


(4.6)  E  Z  X.r.1/a|  <  a~2h '  E(K  -K  )2  n,  < 

K  j  j  I  -  n  n+1  n 


which  is  summable  in  n.  Therefore, 


r_n  x.rT1//a -e i'n  x.(r.1/aAi) 

1  j  j  i  j  j 


converges  in  probability  in  H.  A  short  calculation  gives 


n(logg  n) 


2/a  ’ 


,v.'v  ■  V'.-  V  v  V  V  V  V  ■.>  V 


",  •*.  -r.  \  -V.  •*  -\ 

v*  .%v  -  ' 


\*\*\'*\*V  1 


,,  ...  .  .  ...  A  _  -1/a^.l/a  _  n  ,-,-1/a  .  .  .  f 

(4.7)  A  (a,o)  =  E  a  6  (i  .  A  1)  -  o 

n  1  J  J 


-1/a.l/a  -1/a  .  .  2 
a  o  a  1+r  r 
n 


U 

.  o-1V/0‘f‘n  «t'1/0Al)  -■  C 


-1/a 


l+a”2/a52/a  t‘2/a 


)  dt 


-*■  finite  limit  =  A(a,o)  as  n  -*■  «  . 


Therefore  X^  converges  in  probability. 

Theorem  4.8.  If  (4.3)  is  the  log  characteristic  function  of  a  stable 
law  on  H  then  the  series 


i-l/a51/azr.  (x> r71/a-(EX. )  f j 
1  3  1  1 


-1/a 


"j-1 


Hct"2/a62/Ct  t 


2 /a 


)  dt) 


converges  a.s.  in  K  to  a  random  vector  with  log  characteristic 
function  (4.3). 

Remark.  Centering  is  not  needed  for  the  case  a  <  1,  nor  is  it 
needed  for  the  symmetric  case  which  will  appear  in  [4]. 

Proof.  Since  converges  in  probability  and  X^  is  an  inde¬ 

pendent  series,  we  conclude  by  ([2],  Theorem  5.3(6))  that  X^ 
converges  a.s.  in  H.  Recall  that  with  probability  one  finite 

M  such  that  (n  _>  >1)  =*>  (j  smallest  N’(n)  with  n  =  K^^).  Then 
Vn  >  m, 


3. 


The  first  term  equals  the  first  term  above.  For  the  second  term 


T3  H 


(A. 9) 


v(>Kn))  -l/a.l/a-n  ,-1/a 
X  =  a  a  iL.  . 

1  1  J 


■(EX)  (t  1 )  d  t  : 


(EX  ) A  (a,s) 
1  n 


s““  for  "  -  “•  '  K»(„),">  (aXU>-l  '  r„  i  aS(„>) 


(4.10) 


f  N(n)  , 

I  (t  i/:iAl)dc;  <  (N(n))'1 


By  the  law  of  the  iterated  logarithm,  a.s.  eventually  as  n  -*•  00 


(4.11) 


f L  n  .  —  1  -'ci  r'n  -1/a 

:  (t  a  l)dt  =  t  dt  < 


2  -/n  log  log  n  (n  +  o(n) ) 


.  -n  (n  ■'  r  ) 
n  n 


Therefore, 


(4.12)  X(N(n))  =  d-1  /CX61/0l{ z”  X.:.1/a  -(EX  )  ’  a  i/jla  l ) d t }  - 

1  2  J  1 


A(a,:)  (EX1)  +o(l) 


which  converges  a.s.  in  H  to  a  random  vector  with  log  characteristic 
function  (4.3).  q 

Remark.  lhe  centerings  used  above  also  have  an  interpretation  involving 


(EX^)  E(T  "  1),  which  will  not  be  given  here. 


/  / 


5.  Symmetric  Case,  Multiple  Representations. 


In  this  section  we  do  not  assute  that  is  distributed  according 

to  the  measure  y  ,  or  even  restrict  its  distribution  to  S.  Suppose 
X,  e,  T  are  mutually  independent  sequences,  with  e  ,  "  as  in  section 

3,  and  (X_.,  j  >_  1}  i.i.d.  in  H  with  EljX^Jj  2  < 

Remark.  Series  of  form  I1?  £.X  7  will  be  termed  symmetric. 

-  1  1  1  1 

Lemma  5.1.  The  symmetric  series  2^  e  X  F.^a  converges  a.s.  in 
-  1  J  J  3 

H  and  the  log  characteristic  function  of  its  limit  is  e!(x,X  );^E(u), 

A  r  00  dr 

VxeH,  where  B(a)  *  !  a(cos(r)-l)  — —  . 


Proof .  The  arguments  needed  are  similar  to  those  of  theorem  4.8,  but 


easier  in  the  symmetric  case.  For  each  x  c  H, 


(5.2)  loge  Ee 


.  ,  -  'n  -,-1/a. 

‘("-i  y;i  W*  ' 


r  00  ic  (x.X  )  ;  r  ^ 

Zj  -1/a  (e 

J  a  r 

n 


r  ®  ic. r  .  ,  <  , 

=  aE  ;  (e  -1)I(  (x,X  )  ,  <  ra  (x,X  )  — 

1  —  n  1 

^0  r 


=  a  E 


f  00 

f  (cos  (r )  -1)  I  (  (x,X,  )  ,  <  raJ/:t)Ux,L)  ;  3  ~ 
L  r  -  n  i  -*• 


,  (cos(r)-l )  -  — - 
-  1  -Hi 

>  0  r 


From  ([3],  corollary  2.1),  the  limit  in  (5.2)  is  the  characteristic 


V.  .*. v. 


\  *v  V  \ 


[ 

i 


function  of  a  (symmetric)  stable  lav.’  on  H.  Since  the  sums 

-^n+1  .,  t.-1/ol  .  .  .  .  t 

l  C.X.i.  are  independent  and  symmetric,  we  nave  :rom  -  j  , 

K  i  j  1 

11  V 

n  -l/a 

theorem  5.3(1))  that  1,  e.X.T.  converses  a.s.  in  H.  Since 

1  J  J  J 

n  „  - 1  / ; 

eventually  a.s.  K  , ,  <  K  +1  as  n  -»  00  ,  we  conclude  I,  e.X.7. 

J  n+1  -  n  1  J  J  J 

converges  a.s.  in  H.  q 

Several  series  may  represent  the  same  stable  law. 

*•  ||  oc 

Theorem  5.2.  If  Ej,  X^|,  <  “  then  for  every  x  £  H, 

E|  (x,HL)  !a  =  Ej|  [aE,  (x,X^)  |a  where  X^  is  distributed  on  S 
according  to  the  measure; 

x  >•  ;;a 

P  (X*  e  A)  -  El  (r~ -  C  A)  --  --1  1  ■ 

i  ii  Aii.  r  y  :  a 

1  L'‘l" 


Proof .  For  every  x  C  H, 


X, 


E;(x,X  )'r  =  E.(x,  ™)  *  X 
1  I 


=  Ej^r  EiCx.X^ 


As  an  example  of  the  above,  every  symmetric  stable  law  on  H  has 

°0  —  l/r*  ,  1/r, 

a  construction  Z^.  X^  7  ^  /E(|Z^  )  in  terms  of  an  independent 

sequence  Z  of  i.i.d.  standard  normal  r.v..  Conditional  on  the  seque 
T  and  X  the  symmetric  stable  is  Gaussian.  That  is,  symmetric  stab! 
laws  are  particular  mixtures  of  Gaussian  laws  with  rero  means  and  d if f 
covariance  kernels.  The  latter  will  not  in  general  differ  only  by  sea 
though  this  is  necessarily  true  for  H  =  . 
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6.  Symmetric  Stable  Random  Set  -unctions. 

Let  ■’ t  .  ,  j  _>  1;  be  i.i.d.  taking  values  in  a  measurable  space 
with  measurable  sets  generically  denoted  A.  Suppose  c.  ,  7 

are  mutually  independent,  where  the  latter  three  sequences  are  as  in 
section  5.  Define 


(6.1) 


'  00  -1/u 

x(a)  =  i (t .  ■  A)x.c.;. ,  va 

1  J  J  1  J 


Theorem  6.2.  The  series  (6.1)  is  a.s.  convergent  for  each  A,  is 

jointly  symmetric  stable  for  finitely  many  A  at  a  time,  and 

X(A^)  ,...,  XfA^)  are  f°-  each  n  _>  1  mutually  independent  if 

At,..., A  are  r.utuallv  disjoint. 

1  n 

Proof  .  For  each  n  j  1,  x €  H,  real  numbers  r„,...,r  ,  and  measur 
able  sets  A. , • • . , An , 


i(x,r;  r,  M ( A.  ) )  i  ,  :n  ,  r ,_I ("  c A.  )  ‘  (x , X 

/  s  ~>  \  1  r-  IK  K  _  J  -  1  t.-  1  rv  J  * 

(6.3)  log^  ^  c  “  loge  ^  e 


(5.1) 


=  3  ( u )  £ '  ( x ,  X  )  ~e  :1‘  r,  I  (t  t  A,  ) 


If 


,A  are  mutually  disjoint  then 
n 


(6.1) 


..n 

"  k=  1 


rkI(’l  c  V  '  "k*l 


r  ( " ,  -  a, )  . 


□ 
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Remar  k. . 


The  s  im 


way  ir.  which  a-d  eper.d  er.,  e  ,  d  inei-.s  *  .  r.ai  a:ru  fare  ,  a 
dependence  are  ider.t  if  red  with  r.uf  allv  i:id  crer.d  er.t 

--1/a  ..  . 

sequences  .  ,  X ,  .  . 


Remark .  Schiider  [  8j  and  Ruelbs  [2]  have  explored  a  ret r eser.tat to: 
multidimensional  symmetric  stable  r.v.  by  means  of  a  stochastic  in 
with  respect  to  a  on.e-d mens ional  stable  independent  increments  j.r 
Theorem  (4.8)  and  lemma  (5.1)  sharpen  and  extend  sue:,  rerresent.it  i; 
by  connecting  them  with  the  Fergu son-Kla ss  representation,  making 

explicit  the  choice  of  coefficients  required  to  obtain  each  stable 

4 . 

law  ,  and  establishing  H  convergence  of  the  indicated  series. 

°C  ~  1  / 

Remark.  Suppose  K  •  1  and  V.  =  ■  T,  . '  \  1  •  k  •  urc 

K  1  inJ  KJ  “  — 

independently  constructed  (as  per  (5.1))  symmetric  st  i'-T.  v  r.v.  tan 

values  in  H.  Then  for  an  arbitrary  choice  of  real  numbers  r  ... 

_  ao  ^  ^  —  1  J  ri 

the  sum  ^  r, Y  is  representable  c.T.  '  ‘  where  ■  '  ,  i  ■  1 

1  k  k  r  1  J  J  J  — 

,  „  p  i  '  , 

are  i.i.d.  and  =  an  equ iprobable  random  selection  :rom  K  i 

.1/a 

•  •  •  > K  rn’jni'  uses  the  property  (discussed  in  section  3)  of 

Poisson  processes  run  simultaneous!'.'. 


Including  the  nc: -svnmet  r  ic  stable  laws. 


Stable  Random  Sanctions. 


micable  Star  lunarv  Svm.net  r 


basically,  we  seek 
stationary  random  fund 
char ac t er ist ic  function 


to  construct  the  stable  analogues  of 
ions  having  a  harnori:  decomposition. 


of  such  a  Gaussian  random  function 
i(A_.t.  ) 

rn 

y  re 
k=l  k 


i('\,t.)  ; 

1  K.  O 


Gaussian 

The 

Involves 


~here  t  is  generic  for  a  point  or  the  domain,  and  (*  ,  )  is 

1  o 

a  random  linear  function  on  the  domain.  The  stable  analogues  of 

these  Gaussian  random  functions  have  characteristic  functions  that 
th 

employ  an  a  -power  ir.  this  integral  instead  of  the  2,  but  are 
otherwise  the  same.  Define  Vt , 


(7.1) 


X(t)  = 


:os((' 


.  .  t )  - 
J  ° 


i .  )  X .  c  . 
J  J  2 


.-1/a 


where  '  are  i.i.d.,  2  are  i.i.d.  uniforms  on  [-2,2],  tt.X.T: 
are  as  in  section  6,  and  A,  2,  t ,  X,  7  are  mutually  independent 
sequences.  The  series  (7.1)  is  a.s.  convergent  in  H  for  each  t 
by  lemma  5.1.  The  random  function  X(  )  is  clearly  stationary 
because  3  are  uniform  on  [-0,7],  but  this  will  also  be  a  simple 
consequence  of  the  form  of  the  characteristic  function  which  we  now 
compu  t  e . 


(7.2) 


log  E  e 
e 


1(X’:k=l  rkX(tK)} 


i  7  ,{1?  .  r  cos((A.,t  )  +G  . ) }  c  .  ( x ,  X 

=  log  Ee  ^  k=1  k  J  ^  2  3 


(5.1) 


B(a)E|(x,>1)  I3  E  :"=1  rR  cos(  (/^ ,  t^  +  9^ 


'V  V  W  ■C'JV 


The  final  tem  in  the  right  side  of  (7.2)  reduces  as  follows,  with 


^  rn  1  k  o 

z  =  I  r._  e 


(7-3)  E  i  E1  rk  cos((Alftk)o  +  3^ 


„-ai  101  -  ‘ 

E2  ; ze  +  ze 


,  -Lw/  - 

!  J  QL  — C1„Z  •  ,  -1-  1  | 

*  E  j  z  |  2  E  ,1+z  ze  , 


=  E ! Z ! a  2  2C(3) 


where 


C(a)  =  ; 


in  a  ,  i  ’  Vt  —2  ir  & 

1  +  e  :  dr  =  .  i  +  e  v  x  1  d' 


for  all  real  '!  .  We  have  therefore  proved, 

Theorem  7.4.  The  random  function  defined  by  (7.1)  converges  a.s.  in 
H  for  each  t,  and  has  log  characteristic  function 


i(A. ,t, ) 

,-a  ,  ,  ,  .  irii  1  k  o  ■  a 

2  3(^)C(a)E,I  r.  e 

k-  X  r* 


for  all  n  >  1,  r_,...,r  ,  t  ,...,t  .  m 

—  i  n  l  n 

Corollary  7.5.  The  random  function  (7.1)  is  non-ergodic  for  each  n  <  2 
Proof .  By  using  (7.2)  the  construction  (7.1)  remains  valid  if  t.  are 


standard  normal  sequence.  Conditional  on  the  sequences 


1  <  2. 

are 

& 

i. i .d  . 

< 

T,  the 

1 

>: 

V 

5 
.*  \* 

- \  %s,  ■ 

-  n- 
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8 .  Operator  Stable  and  Infinitely  Divisible  Multidimensional  Laws. 

Infinitely  divisible  laws,  of  which  the  operator  stable  laws  are  a 

special  case  with  particularly  interesting  structure,  are  treated  in 

Part  II.  In  brief,  this  is  what  happens:  A  construction  of  infinitely 

divisible  random  vectors  is  given  by  ,X^  )-y  },  in  which  the  real 

function  H  is  monotone  decreasing  and  positive  for  each  value  X.,  and  is  deter 

mined  from  the  Levy  measure.  A  construction  of  full  operator  stable  random 

vectors  in  a  finite  dimensional  real  vector  space  is  I{A(r.^)X  -Y  },  in 

3  3  3 

which  the  vectors  y  are  non-stochastic  centerings,  { A ( t )  =  exp(Blog  t) ,t  >  0} 
is  the  group  of  linear  transformations  figuring  in  the  definition  of  operator 
stability  (e.g.  Sharpe  IgJ),  and  the  vectors  X  are  i.i.d.  from  a  probability 
measure  (a  factor  of  the  Levy’  measure)  on  a  set  of  generators  of  the  subgroups 
induced  by  A.  The  methods  of  sections  4  and  5  carry  over,  as  will  now  be  in- 
dicated.  If  X  is  any  i.i.d.  sequence  in  R  ,  and  X  is  independent  of  7,  then 
VxeRd,  n  >  1, 


(8.1)  log 


E  -  f" 


E(ei(x’A(t)Xl}-  1)  . 

a  t 

n 


As  usual,  the  symmetric  case  is  simplest.  If  we  examine  Sharpe's  Theorem  5, 
we  discover  that  the  limit  of  (8.1)  is  precisely  the  form  taken  by  the  oper¬ 
ator  stable  in  this  case,  provided  we  choose  for  the  distribution  of  the 
probability  measure  figuring  in  Sharoe's  representation  of  the  Levy  measure 

as  a  mixture,  this  measure  being  placed  on  (Sharpe's  notation)  generators 

B  _  1 

0  characterized  bv  s  Mg',  t  ?  :  t  >  s  /  =  s  ,  V  s  >  0.  Arguing  as  in  section 
00  1  d 

5,  we  conclude  7^  A(7  ^-j  converSes  a.s.  in  R  and  has  the  log-charac teris- 

tic  function  which  is  the  limit  of  (8.1). 


■  V  V  V  V  ’/  *.  V  V  W  V  '  *  *.*  V  V  v 


*»  »  »■  *  •»  » 

« • . 

»  ^  •  ..  _  *  _  « 
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9.  Priority  of  P.  Lew 


P.  Levy  has  anticipated  the  series  constructions  of  one  dimensional 
stable  r.v.  with  a  <  2.  For  the  case  of  a  positive  stable  with  a  <  1, 


to  scale  and  location,  this  construction  is 


is  £  F.  ,  with 

.  j=1  J  .  +. 


{I\,  j  >  1}  being  the  arrival  times  of  a  Poisson  process  (on  R+)  having 

unit  intensity  function.  Levy  writes  the  series  in  the  form  T  U  ,  where 

L  x 
X 


(9.1) 


(9.2) 


(U  ,  x  >  0}  are  independent  r.v.  and 

X 


P(Ux=>:>  -  ^^x"  0) 

X 


Here  is  my  abstract  of  the  key  parts  of  Levy's  (1935)  arguments  for  the 
above  case: 


(9.5)  [  =  x  a  <“]=*>  [{u  :  U  4  0,  x  >  x  }  is  finite  for  xn  >  0 

i-ra  u  x  x  u  0 

JX  0  X 


and  also. 


(9.4)  [ 


0  +  o) 

0  1-“  0 


IE  l  V  i  0] (as  x„  +  0) 

-  x  U 

x<  x0 


Therefore,  for  arbitrary  c^,  c0  >  0  (defining  c^  =  c“  +  c^  and  taking 
independent  copies) , 


(9.5) 


:.  I  U(1)  +  c.  I  U( 
1  L  x  2  L  x 


(2)  Dist*  „  (i)  „  (2) 

y  uw  =  ,  \L  +  >y 

L  X  X  L  X 

X  X 


Dist . 


Dist . 


;  iM 
x 
x 


c,  /  U  (=*■  stable) , 
3  L  x 


where  (Yv  ,  x  >  0}  have  respective  intensities  c~  a  dx/x  ,  k=  1,2,3, 

X  K. 

and  are  independent  for  k  =  1,2. 

The  above  arguments  do  yield  a  proof  of  the  representation  if  we 

apply  them  to  the  independent  sub-sums  £{U  :  x  G  b^]},  n  >_  1, 

-a 

where  b  =  log  n.  This  is  essentially  the  argument  of  Ferguson-Klass 

(1972).  The  particular  choice  of  b^,  n  1  is  one  which  ensures  that 

eventually  as  n  -*•  00  each  sub-sum  contains  at  most  one  summand,  so  it 

really  is  (almost)  as  though  one  could  add  independent  U  one  at  a  time 

toward  x  0.  A  quite  different  justification  is  to  interpret  7  l* 

x 

as  a  generalized  process  driven  by  "white  noise"  {U  ,  x  >  0} . 

X 

Levy's  observations  are  easily  overlooked.  Ferguson-Klass,  Vervaat 
(1979),  LePage-Woodroof e-Zinn  (1979)  (in  manuscript  form),  rediscover 
the  Levy  construction  as  byproducts  of  the  following  independent  pursuits 
respectively:  (F-K.)  -  representing  the  positive  non-Gaussian  part  of  an 

independent  increments  random  function  as  the  sum  of  its  ordered  jumps. 

(V) -examining  a  shot-noise  associated  with  the  asymptotic  behavior  of  the 
solution  of  a  stochastic  difference  equation  as  time  is  increased, 

(L-W-Z) -studying  the  limit  behavior  of  the  normalized  order  statistics 


* 


v  -IS- V V  •  -  .  .  .  .  .  .  .  .  --  . 


from  a  distribution  attracted  to  a  stable.  Resnick  (1976)  reconciles  the 
Ferguson-Klass  construction  with  the  Ito  representation,  meaning  by  the 
latter  Ito's  generalization  of  Levy's  stochastic  integral  construction 
by  a  Poisson  random  measure. 

Acknowledgement .  The  estimate  (4.6)  was  suggested  by  J.  Zinn.  Helpful 
discussions  were  also  held  with  V.  Mandrekar,  M.  Woodroofe,  and  M.  Steele. 


I 

I 

References 

[1]  Ferguson,  T.  and  Klass,  M.  (1972).  A  representation  of 

independent  increments  processes  without  Gaussian 
components.  Ann.  Math.  Statist.,  43,  5,  1634-1643. 

[2]  Hof  fman-J«5rgensen,  J.  (1977).  Probability  in  B-spaces. 

Lecture  note  series.  No.  45,  Aarhus  Universitet. 

[3]  Kuelbs,  J.  (1973).  A  representation  theorem  for  symmetric 

stable  processes  and  stable  measures  on  H.  Z.W.  Verw. 

Geb.  26,  259-271. 

[4]  LePage,  R.,  Woodroofe,  M.,  and  Zinn,  J.  (1980).  Convergence 

to  a  stable  limit  distribution  via  order  statistics. 

(5J  Levy,  P.  (1935).  Proprietes  asymptotiques  des  sommes  de  variables 
aleatoires  independantes  ou  enchainees.  Journ.  de  Math., 
tome  XIV.-Fasc.  IV. 

[6]  Logan,  B.F.,  Mallows,  C.F.,  Rice,  S.O.,  and  Shepp,  L.A.  (1973). 

Limit  distributions  of  self-normalized  sums,  Ann.  Prob., 

1,  5,  788-809. 

|7]  Resnick,  S.  (1976).  An  extremal  decomposition  of  a  process  with 
stationary,  independent  increments?  Tech.  Rep.  No.  79,  Dept, 
of  Statistics,  Stanford  University,  Stanford,  California. 

]8]  Schilder,  M.  (1970).  Some  structure  theorems  for  symmetric  stable 
law.  Ann.  Math.  Statist.,  41,  412-421. 

19]  Sharpe,  M.  (1969).  Operator-stable  probability  distributions  on 
vector  groups.  Transactions  Amer.  Math.  Soc . ,  135 ,  51-65. 

J10]  Vervaat,  W,  (1979).  On  a  stochastic  difference  equation  and  a 
representation  of  non-negative  infinitely  divisible  random 
variables.  Adv.  Appl.  Prob.,  11,  750-783. 


4ft  l  J 


»Vi y.iiW*. 


V  1 


■  .  ■ 


1* 


23 


